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sound speed, we flnd a lower limit on the sound speed related to the amplitude of the 
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1. Introduction 

Inflationary model building has been a popular pursuit in the years since data from the 
cosmic microwave background became precise enough to begin distinguishing scenarios. 
Model builders have taken either the approach of attempting to constrain the terms in 
a "generic" model from observation ('reconstruction') or have searched for particular 
models with interesting features in the context of supers jmametry, string theory or 
related ideas. However, it is still not clear what range of observables, particularly for the 
amplitude of tensor fluctuations and non-Gaussianity, is reasonable or expected from 
fundamental theory. In addition, 'reconstructors' often proceed without clearly stating 
or exploring the assumptions and context for the theories they are reconstructing. 

In either approach, there are very basic guidelines coming simply from the consis- 
tency of the perturbative approach to calculating fluctuations. The implicit assumption 
is that the inflaton obeys an eflective fleld theory (often leaving aside the important 
question of vacuum selection, which we too will ignore here) and on scales below some 
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cutoff its ffuctuations can be calculated order by order. The conditions for this per- 
turbative expansion to be valid are trivial and well known in the context of slow-roll. 
However, smooth, single-field slow-roll inflation cannot produce large non-Gaussianity 
0] . As an additional complication, single field slow-roll models are surprisingly hard 
to find in string theory and do not seem to contain features that shed light on their 
fundamental physics origin |]^, ^. Given all the interest in other models, then, it is 
interesting to explore the extension of the consistency constraints to other cases. Here 
we do so for models with small and changing sound speed like k- inflation ghost 
inflation and Dirac-Born-Infeld (DBI) brane inflation [0. We also look at a toy 
example with multiple fields. This is a simple idea, and some of our results have been 
obtained in some form previously [§, |^ so we emphasize the generality of the approach 
and some interesting applications. 

In the standard picture of inflation, the classical behavior of a scalar field drives 
the background accelerating expansion of the universe and quantum fluctuations of 
the field during inflation are eventually imprinted as perturbations in the gravitational 
potential after inflation. There are two basic consistency conditions here: first, that the 
energy density of fluctuations is not so large that it overwhelms the background, and 
second that interactions of the field are small enough that the amplitude of fluctuations 
can be well calculated from the term in the action quadratic in fluctuations. We will 
call the first condition the gradient energy condition, and in slow roll it is satisfied if 
H/Mp < 1 which is anyway a good condition that keeps the energy density driving 
inflation sub-Planckian. ^ The second condition (the interaction picture constraint) is 
never violated in slow-roll as long as H/Mp < 1 and the slow-roll parameter e is less 
than one (also the condition for accelerated expansion). 

Particularly as we become interested in models that have large interactions gener- 
ating large non-Gaussianity, it is useful to put these constraints in the slightly more 
formal, but also most general, language of the action. Schematically, the action for a 
general scenario can be written in powers of the fluctuations 

^ = S-o + ^2 + ^3 + • • • (1.1) 

where 5*0 describes the classical evolution. Si is absent when the background fields 
satisfy their equations of motion, 5*2 is the quadratic piece that gives the two-point 
function (the power spectrum), and higher orders are interaction terms. Then the two 

^The gauge invariant method to calculate backreaction effects was discussed in and we tliank; 
the referee for pointing this out. Note that the condition we give in this paper does not include the 
stochastic behavior of the field and therefore is strongest near the onset of inflation. In that sense it 
is a minimal condition. 
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conditions above can be phrased schematically 



So > S2 (1.2) 

S2 > Sn, n > 3 

where the amplitude of terms can be estimated using the the average amplitude of 
fluctuations near horizon crossing. Since calculating higher orders in perturbation 
theory is a long chore, we will focus on the first case of the second line above, 5*2 > 5*3. 
This is really a statement about a loop contribution to the two-point function, where 
we have suppressed a log factor. 

The main application of these ideas in this paper is a lower bound on the sound 
speed for which the perturbative expansion is valid. This bound is particularly in- 
teresting since it imposes a strong constraint on models like DBI inflation, where the 
sound speed can vary considerably. Given the sound speed at some scale /cq, we can 
compute the number of e-folds obtained in the perturbative regime, as a function of 
the running. This also has implications for bounding the expected non-Gaussianity on 
non-CMB scales, where there is as yet no interesting observational constraint. This 
bound coincides with the expression found in by an alternate method. This bound 
also implies that eternal inflation is outside of the pertubative regime in any model with 
small sound speed, which was realized independently in a recent paper of Creminelli et 
al [0. Finally we also find a bound on the number of fields, A^, in a simplistic multiple 
field example which agrees completely with recent bounds from renormalization of the 
Planck mass fill and from black hole physics [l^ applied in de Sitter space [[1^ . 



In the next section we review the basic framework of the perturbative calculation. 
In Section |^ we examine in detail the regime of validity of the calculation in the case of 
slow-roll and general sound speed models with scale-dependence. Implications of the 
bounds for scale-dependent models are discussed in Section |[ Multi-field inflation is 
discussed in Section ^ and Section ^ concludes. 



2. Set-up, Gauges and Conventions 

Here we outhne the class of theories we consider and choose a gauge for the calculation. 
We will first consider the following system 

S=IJ d'x^g {MlR + 2P(X, 0)) (2.1) 

where is the inflaton and X = —^g^'^d^(f)du(p- In usual slow-roll, P{X, 0) = X — V{(j)). 
We are interested in an inflationary background where the inflaton is homogeneous with 
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perturbations t) = (j)o{t) + 6(j){x, t). The background metric is the usual Friedmann- 
Robertson- Walker (FRW) metric 



ds'^ = -dt^ + a^{t)5ijdx'dx^ (2.2) 

Einstein's equation together with the continuity equation completely describe the (in- 
flating) background 

p=-3{p+P) (2.4) 

where p is the energy density, P is the pressure and if = ^ is the Hubble parameter. 
It is useful to define the speed of sound: 

dp P^ + 2XP^x 



and the "slow-roll" parameters 



H XP 



X 



if2 M^H"^ 



e 

7h 



(2.6) 



where Px means the derivative of P with respect to X. To analyze the validity of the 
perturbative analysis, we will work (following 0) in the ADM formalism, where the 
metric is Jl^: 

rfs^ = -N^dt'^ + hij{dx' + N'dt) {dx^ + N^dt) (2.7) 

The lapse functions, A^, N'^ function as Lagrange multipliers which can be eliminated 
through their equations of motion. Since the metric also fiuctuates, it is important to 
choose a gauge and correctly account for all the fiuctuations. In this paper, we will be 
using two different gauges, and we write only the scalar part: 



1. In the CO moving gauge, 50 = and hij = a?e^^5ij. One can solve for the lapse 
functions in this gauge P to find to first order A^* = di {^o- '^jj + d'"^ ( 2H^cl ) ) 

and N = 1 + jj. is constant outside of the horizon and corresponds to the 
measured curvature perturbation in the post-infiationary era. 
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2. In the spatially flat gauge we set hij = a^5ij and expand the action in terms of 
5(j). The lapse functions are then N = l + ^P^xH and A^* = did-^^f^ (^-^6 



16|. In this gauge the amplitude of fluctuations of the inflaton is not constant 



after horizon exit and will change. 

One can move from one gauge to the other by a gauge transformation which to 
first order in the fluctuations is 

C = -f ^0 (2.8) 

Working in the comoving gauge, we can expand the action in fluctuations of the metric 
(. The first order vanishes since it is proportional to the background field's equation 
of motion while at second order 

^2 = J dtd^x (^a^-^C^ - a<dCf^ (2.9) 
Terms proportional to three powers of ( {S3) were calculated explicitly in P] for usual 



slow-roll and for general sound speed in ||T7| , p!8| , p!5|] . The full expression is quite long 
but contains the following significant terms: 

^3 = [ dt£x^{e -2s + 1- c^)C(9C)' + . . . (2.10) 

The last three terms in the parentheses cancel for a canonical kinetic term where = 1 
and s = 0, and the first term is the usual slow- roll piece which carries factors of e^. 

2.1 The Variance of C, 

While the system is stable at the classical level, quantum perturbations will give a 
non-zero variance to the random C, perturbations. In calculating this variance, there 
are two assumptions involved: 

1. The FRW background is valid and the correction from the perturbations is small 
compared to the background. 

2. The dynamics of the C, field is dominated by its quadratic action and higher order 
interactions are small. 

Therefore, as usual, one solves for the perturbation C, using 82- The perturbations are 
decomposed into momentum modes 

uik,t) = [ d:^xCix,t)e-'^-^ (2.11) 
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and ( is quantized 

Cik,t) =u(k,t)a(k)+u*{-k,t)a\-k) (2.12) 

where as usual [a{k), a\k')] = (27r)^5^(/c — k'). The solution to the equation of motion, 
assuming all three slow-roll parameters above are small (e, 77, s -C 1), is 

Uk = (1 + ikc,T)e-'^'^^ (2.13) 

where dt = adr is the conformal time (t ~ — ^ for e ^ 1). Before horizon exit 
Cgk > aH and oscillates with a decreasing amplitude (oc r). After horizon exit \uk\ 
is constant and on the order of 

Uk ^ (2.14) 

The two-point function is defined 

|C(fci)C(^2)) = {2Tif5\k^ + k2)Pc (2.15) 

= {2Tlf5^{ki + k2)27T^Vck~'^ 

where is the dimensionless power spectrum (variance), and assuming the spectral 
index doesn't run "P^ = A{kQ){k/kQ)""'~^. Using the solution above, we also have 

(C(^i)C(4)) = i2n)'6'ih + hMk)]' (2.16) 
i2nf6%h + k: ^ 



Averaging over a sound horizon volume, C{x,t) can be thought as a random variable 
with mean zero and fluctuations of size „ ^ — . Note that a smaller sound speed 
means a smaller horizon for scalar modes. From the variance of ( we can do a gauge 
transformation and compute the variance of S(f), we find 

(2.17) 



where we have used the definition of e in Eq.( p.6|) and 2X = 0^. For a canonical kinetic 
term the variance is of order H'^/An'^. Now let us check under what conditions our two 
basic assumptions are valid. For the rest of this paper, we will be mostly interested in 
order of magnitude calculations and we will drop the factors of 2n and take ( ~ 
over scales of length L and time of order 1/H. 
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3. Perturbative Limits During Inflation 



The first condition one must have for a valid perturbative analysis of inflation is that the 
background is stable. That is, we demand that the contribution to the energy density 
coming from the pertubations is smaller than the smooth constant energy density that 
drives inflation. This is ensured by demanding that 5*2 < 5*0. For an inflating solution 
where the potential dominates, one gets that 5*0 ~ V{(f)(_o)) which is of order H'^M^ 
using the Friedmann equation. Given that ( ~ over a Hubble scale with 

— r^H A _ ^ (3 1) 

dt ' dx Cs ^ 

the gradient energy bound for a general single field model is 

S2 ea'\dCf 

So ^ Hm^ ^ ■ ' 

< 1 (3.3) 



M2c3 

•p s 

The kinetic energy term in 5*2 contributes at the same order. Note that for usual slow- 
roll models with Cg = 1, this condition just reduces to having a subplanckian Hubble 
scale. When Cg is small the gradient energy induced from the perturbations is enhanced 
in two ways: the variance increases while the horizon decreases. We therefore find that 
H must be significantly less than the Planck scale in a small sound speed model. 

p 

This condition can also be stated in terms of the stress-energy tensor, using the spatially 
flat gauge. The gradient energy condition is 

^ < 1 (3.5) 

The first non-zero contribution to (5Tq ) contains the same terms (kinetic energy and 
gradient energy of fluctuations) as the action at second order in fluctuations, 5*2 given in 
Eg. ( p. 91) . For potential energy dominated inflation, we can shorthand Tq as SM^H^ ^ 
V{(j)), which is also the dominant term in 5*0. 

3.1 The Expansion Parameter to Third Order 

While the first condition is somewhat trivial, the condition demanding that all the 
higher order interaction terms are small is more interesting. Higher order terms in the 



-7- 



action will correct the two point function. An example loop calculation was done in 
(for a case with multiple fields, which we return to in Section ^ with the conclusion 
that, up to a term logarithmic in scale k, ^ 1 is required for the loop contribution 



to be small. The importance of the log term is considered in detail in |2T], |23 
but we emphasize that the simplest conditions must hold when the log is of order 1. 
Although for a generic model the loop calculation is quite involved, the bound that 
results is not surprising. The well-studied local model, for example, parametrizes non- 
Gaussianity in a simple way that helps us see what to expect. In the local model the 
primordial curvature perturbation in real space is written as a Gaussian piece, (g, plus 
a quadratic correction 

C = C, + lfNL{Cl-{Q) (3.6) 

A quick computation shows that the two-point function is not significantly corrected 
while f^j^ (Cg) <^ 1 and in that case the fluctuations are nearly Gaussian. For fluc- 
tuations generated during inflation, 'nearly Gaussian' means that the primary infiaton 
field is nearly free. In single field slow-roll this statement is related to the flatness of 
the potential. More generally, 'nearly Gaussian' means that we can solve the quadratic 
equation of motion for the fluctuations and treat interactions as perturbations, that is^ 

Sn/S2 « 1. 

With these examples in mind, we assume that the requirement S3/S2 <^ 1 also 
follows from the loop calculations with the terms we are interested in. Requiring 
5*3 < 5*2 and using the 1-sigma size of ( gives 

^r^^{e-2s + l-cl)<l 



clMpiec.y/' 



(e-2s + l-cf)<l (3.7) 



Depending on whether the model in question is a slow-roll case (where Cg = 1 and 



quick, approximate evaluation of the correlation functions can be made assuming the dominant 
contribution comes from one Hubble time near horizon exit. Then (C") « At{Ax)^ C2{Cn / ^2)1^^^^ 



where £„ is the Lagrangian at nth order in fluctuations. Using Eq.(p.9[) and Eq.(3.1), the product 



of the first three terms is about 1, so the dimensionless quantity /T?"/^ ^ £„/£2, which is the 
same as Sn/ S2 as we are using it. The same dimensionless combination of the n-point function scaled 
by the two-point appears in the expansion when one writes the probability distribution function as a 
Gaussian modified by a series of cumulants. 
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0) or is a small sound speed case {1 — cj. > e — 2s) there are two different bounds: 

< 1 slow-roll (3.8) 



M2 

p 



Jj2 

< cl small Cs (3.9) 



When the sound speed is small and Cs < 1/e (always true during inffation), this bound 
is more constraining than the gradient energy bound. This last equation is particularly 
interesting since one can write it in terms of the observed magnitude of the curvature 
power spectrum to find a lower bound on the sound speed 

This gives a direct experimental bound on how small Cg was when CMB scales left the 
horizon. Indeed, Cosmic Microwave Background observations indicate "P^ ~ 10"^ ||32|| , 
implying Cg > 10~^/^ for a perturbative analysis of inflation to make sense at CMB 
scales. This relation agrees with the bound found in ||^ by a different method. 

The bound on Cg translates into a (somewhat model dependent) bound on the 
magnitude of non-Gaussianity coming from the three-point function: Iftl/] ~ ^ < 
10^/^. Here |/^^| is calculated at the equilateral limit {ki = k2 = k^) of the momentum- 
space three-point function. The bound again agrees with our expectations from the 
local model. This is much weaker than current experimental constraints at CMB scales, 
which give roughly |/^^| < 300 for small sound speed models However, these 

bounds must be satisfied over the full duration of inflation in order to be able to trust 
the perturbative analysis. For models where the sound speed decreases during inflation, 
we will use this result to derive a constraint on Cg in terms of the required duration 
of the inflationary phase. In particular DBI brane inflation has a sound speed that 
decreases during inflation as the brane moves along a space with continually decreasing 
warp factor. 

We have concentrated on some particular terms in 5*3 (largest when Cg is small), 
but there are certainly others. In fact, there are two families of terms that introduce 
new parameters at each order in the expansion of the action. The first family comes 
from interaction terms in the potential and can be arranged in terms of the Hubble 
slow-roll hierarchy parameters (always assuming the potential energy dominates) [p5| , 
where new parameters enter that are related to consecutive derivatives of e. The second 
is an analogous family coming from the kinetic terms. For example, at third order it is 
useful to define 



X — X^P XX + -X^P,xxx 



6 \ci / 3cg 



(3.11) 



- 9 - 



From the reconstruction point of view, there is a constraint on the new parameters at 
each order. For example, the term in 5*3 containing A is of the order '^}I^3/2 and this 
must also be less than the contribution from 82'- 



.3/2 



< HMp (3.12) 



Note that throughout this section, we have assumed there is no magical cancellation 
among terms at any order. From an effective field theory viewpoint, such cancellations 
would require fine-tuning. 

3.2 Implications for Eternal inflation 

Eternal inflation will occur in any region of space where the magnitude of the density 
perturbation 

^ . ^ 

~ -T7-. — 3.13 

is of order 1, which requires curvature perturbations of order 1 (see for a more precise 
definition). For single field slow-roll inflation, this statement comes from demanding 
that the quantum perturbations of the infiaton dipq ~ H are of the same order as the 



classical motion 50^ ~ This is clear since using Eq. ( p. 81) 



C = -"h, ~ (3.14) 



Now interestingly the bound (|3.1CI|) implies 

< < 1 (3.15) 

and therefore eternal inflation is not perturbatively possible in this class of models. 
It has been claimed in the literature that eternal inflation should actually be more 
common in the phase space of configuration for models with small sound speed. This 
was observed in ^ based on the stochastic approach to inflation and was attributed 



to the increase of the diffusion coefficient thanks to the small sound speed. We have 
shown here that actually all that eternally inflating part of the phase space is outside of 
the perturbative regime and the two-point function cannot be calculated reliably using 
perturbative methods. In the specific case of DBI inflation in string theory, it has been 
verified that for the parameter range allowed by the theory, eternal inflation does not 
occur [^. We emphasize that our statement does not rule out eternal inflation in small 
sound speed model. Instead, we are saying that such an eternally inflating regime will 
always be outside of the perturbative regime. Note that the bound ( p.lOj ) is only vahd 
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for small Cg. In the limit where Cg = 1 the large terms leading to ( p.lO|) vanish as 
(1 — c^) and one is left with the terms that are supressed by an extra power of e giving 

n < ^ (3.16) 

and clearly for small enough e, eternal inflation is possible in the slow-roll Cg = 1 
perturbative regime. While slow-roll eternal inflation can be perturbative locally, there 
might still be regions of the eternally inflating universe where the curvature of space is 
greater than 1/e. Indeed, a global picture of an eternally inflating universe most likely 



requires a non-perturbative description even in the slow-roll picture |29 . 
3.3 Higher Order Terms 

In the previous subsection, we found the conditions for 5*3 < 5*2 and ^2 < Sq. These 
conditions matched our physical expectation in the simplest slow-roll cases, and we 
expect that we have uncovered the 'expansion' parameters. The calculation to fourth 
order has been recently pursued by different groups. To get the action to fourth order 
in the perturbation, one needs to expand the lapse functions N and A^* to second order 



in the perturbations (see for more details on this point). This calculation for the 



action ( p.l|) in the slow- roll limit was carried by Seery, Lidsey and Sloth |^| who found 



that the leading terms (like the gradient terms in S2) are not suppressed by the slow- 
roll parameters and 5*4 scales like H^/Mp. More recently, Arroja and Koyama 



have argued that tensor modes must also be included and give a contribution of the 
same order ^. We see that S4 will be small once the two conditions uncovered in the 
last section are satisfied. While the terms naively order themselves correctly in powers 
of C, a strong scale dependence of the parameters can change that unless a condition 
on the coefficients appearing in the action is satisfied. As we will show, this can be 
particularly troublesome in models with scale dependent sound speed. 

In the small sound speed limit, the leading terms are coming from the self- interaction 
of the inflaton and one can neglect the coupling to gravity [^T], |T^]. Since the slow- roll 
scaling of terms with powers of e is clearest in the spatially flat gauge, where the fluc- 
tuations of the inflaton field, = 0(o) + 50, are non-zero we will work in this gauge in 
this section. We can then expand P{X, 0) in terms of (50. We will not uncover all the 
terms - in particular, we will not display terms proportional to metric fluctuations in 
this gauge, which enter 5*3 and higher. The purpose here is to show the behavior of 
the terms which, from experience so far, give the dominant contribution at each order. 
We will use the result as a starting point to discuss the case of scale-dependent sound 
speed. We switch to Lagrangian notation to avoid carrying around integral signs. 



^We thank F. Arroja and K. Koyama for clarifying this point to us. 
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We examine the kinetic terms, finding at each order new combinations of higher 
derivatives. For reference, two previously defined useful combinations are 



XPx eH'^Ml 



(3.17) 



A — X'^Pxx + -X^P^xxx 



2 

3' 



Notice that X and Px are positive and Cs > 1 would mean P^xx < 0. As discussed in 
33| , there is good reason to believe that this is not consistent with a unitary effective 



field theory. A useful expression is 



1 



2XP 



XX 



P 



X 



(3.18) 



Following |31|, we can Taylor expand the action for the infiaton Pi^X^cj)) in terms of 
the perturbation 50. For brevity, we keep only the Px derivatives, which give a larger 
contribution in models with small and slowly changing sound speed anyway: 



P. XX + ... 



^P,xx 



'-a-2(V50)2] + ^Pxxx0=^5'/ + ... 
o 



-P,xx[^<P -a~\W5<PYY + -Pxxxc. 

4 



- a-\V8(t)f 



(3.19) 



+i^P,xxxx^ 



+ ... 



Inserting powers of 0^ = 2X (the background spatial gradient is negligible), and using 
Eq.( ^.17| ) for the variance of 5(j) and to label the combination gives 



c,P 



Px 1 - ^ + 2XP 



,xx 



+ 



31 

a-- 



2H^ 



3/2 



62c,P 



.X 



X^P 



,xx 



-I 



—p? 

2 



P 



1/2 r 



2d 



X'P 



2\c 



4 

'3' 



+ -X^P 



,xx 



xxxx 



' 3(1 - ci 
2 



ct 



-2a^T.V^ + ... (3.20) 



2 
3' 



-X^P 



XXX 



+ ... 



1-- +4X3p^xx 1-^ 



+ . . . 
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From the last line, we can see that although there will be a new parameter involving 
P,xxxx, the first contribution will scale like C2T^c/'^s where ~ j^ji^^ ~ xPxc ' 
the special case of DBI, the pressure is given by 



(3.21) 



where /(0) depends on the higher dimensional background the brane is moving in. 
Then 



cl = l- 2X/(0) 
1 



P,xx 

P,X" 



(3.22) 



Cs 

m 

(2n-3)!! 



{P,xY-' 



n > 3 



with (2n - 3)!! = (2n - 3)(2n - 5) ■ ■ • 1. Using Eq.fglTl) and Eq.fglsD in the last line 
gives 



X"Px" = (2n-3)!! --1 



1 Sc? 



2("-i) 



,n>2 



(3.23) 



Then rewriting Eq.( p.20| ), we get (using the simplification for the DBI case that 2Ac^/S 



4^3 



Co oc 



~2d 



(3.24) 



V. 



1/2 



-3 ^ -^2-^3 5- 



n-2 



where the Ln are constants. We remind the reader that this only considering the leading 
terms in the action for the small sounds speed and this is specific to DBI in the sense 
that we have assumed a specific form for the parameters A (and equivalent higher order 
generalizations). Nevertheless, we expect the scaling argument of the next section to 
apply in principle to the general case. 
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4. Scale Dependence 



We can apply the bounds from the previous sections in an interesting way for mod- 
els with scale-dependent sound speeds (which are observable as scale-dependent non- 
Gaussianity) . Just as the Hubble parameter evolves slowly in time in a generic inflation 
model, many scenarios may also have a changing sound speed. In analogy to the spec- 
tral index (and using the parameter s defined in Eq.(^.6|)), we can model this by writing 



Csih) { ^ 



(4.1) 



Let us suppose s is a (nearly constant and small) parameter. If s < 0, the sound 
speed will decrease on small scales (and so the magnitude of non-Gaussianity will 
increase). An early discussion of varying sound speed can be found in |Q. To connect 
more directly with observation, we can write the observable magnitude of the non- 
Gaussianity as 

k 



eff 
NL 



kr 



njvG-1 



(4.2) 



where the 'effective' (eff) label indicates that this applies to more than the local model. 
Note also that this parameter f^l is really an approximate amplitude for the full 
bispectrum that is evaluated in the equilateral limit where all momemta are equals. 
In small sound speed models, f^/ oc l/cj. so n^vc ~ 1 = ~2s and there is a model 
dependent coefficient relating f^liko) to c~'^{ko). 

First, we derive a simple relationship between the number of e-folds obtained while 
the perturbative expansion is sensible at this level and apply it to the case of DBI 
inflation. Rewriting Eq. (|3.1CI| ) to display the scale-dependence gives 



Vc<ct 



(4.3) 



A{ko) 



log 



k 
ko 
Ajkp 

ctiko 



< ct{ko) 



kn 



4s 



< (4s -Us + l) log 



k 

ko 



If 4s > 0.05, the sound speed increases with scale and the bound is satisfied for all 
k > kf) (that is, for the duration of inflation). Note that ng — l depends on s in general 
p5| . Otherwise, for H nearly constant, this says that the number of e-folds obtained 
before the bound is violated (or before s or Ug must change) is 

1 



cHko) 



(4.4) 
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0.95 and Aq = 2.22 x 10"9 at 
/cq = 0.002Mpc~^ ||32|, both in terms of initial sound speed Cs(fco) (Eg. ( [4.1|) ) and in 



Contour plots below show this relation, assuming Ug 
-1 



terms of an effective mae; nitude f^/l (Eq.Q), found usine; the dominant contribution 
in DBI evaluated at the equilateral triangle limit. 




-2.0 -1.5 -1.0 -0.5 0.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0^ 
Log[c,(^o)] Log[|/tL(^)|] 

(a) (b) 



Figure 1: (a) Contours of number of e- folds obtained before the bound is violated, as a 
function of the sound speed at some initial scale /cq and the scale-dependence s. (b) Contours 
of number of e-folds obtained in the DBI model before the bound is violated, as a function 
of f^l « — 0.32/c^ at scale ko and with scale-dependence hng- 



The best case scenario for scale-dependent non-Gaussianity to be observable using 
upcoming data on cluster number counts and the galaxy bispectrum has both large 
non-Gaussianity on CMB scales and moderate running {—s ~ 0.1 — 0.3) [^. The 
plots above show that this case is consistent (since the scales differ by about a factor 
of ten, or a few e-folds) as long as s ^ on scales just below the cluster scale, which 
is probably out of the easily observable range. 

In brane inflation where a brane falls down a warped throat, s is roughly constant 
while the brane is in the AdS part of the throat. We can see from panel (b) of Figure 
^ that under these conditions with observable non-Gaussianity at CMB scales and 
running greater than about ung — 1 ~ 0.2, this model cannot achieve enough inflation 
and remain in the perturbative regime. However, there is a possible solution which 
suggests a solution for general sound speed models: the throat is not singular and in 
fact s is zero in the smoothed out tip region where the warp factor is a constant. This 
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feature has been analyzed previously in brane inflation models (for example, [^). In 
order to get enough inflation in this scenario, one can try to find a particular choice of 
model parameters where Cg does not get too small before entering the tip region. We 
can estimate the number of e-folds in the tip region in a simple case, taking H = Acj) 
and using cf) = —2MpH'cs (which can be derived from the continuity equation). Then 

N':^= dct^=-—-J -—d<j, (4.5) 

J^t (p 



where Cs{(f)t) is the (constant) value of in the tip region, 0j is the value of at the 
edge of the tip region and we have assumed (pj (the final value of the infiaton) is much 
smaller than (pf Since (pt < Mp is required by the set-up, we see that the number of 
e-folds at the tip is bounded to be less than iVg*^ < l/4cs(0j) and in general it will 
be much smaller. Then, the number of e-folds from the AdS part of the throat (e.g., 
read off the plot above for a particular case) and the tip region must together provide 
sufficient inflation. This can be imposed more precisely for specific scenarios and may 
be quite restrictive. 

We now turn to a more speculative condition on the coefficients of terms at each 
order. This is speculative in the sense that it assumes the terms calculated in the 
previous section have captured the scaling behavior of all of the most important terms 
that appear when the sound speed is small. In addition, it assumes that one is interested 
in (or expects) a non-Gaussian distribution which can be expanded around the Gaussian 
in a series of converging cumulants. While this is not a necessary feature of non- 
Gaussian distributions, it is commonly assumed by cosmologists and is likely to be true 
when the physics of the fundamental infiaton is simple (e.g., not eternally inflating). 

If the background solution gives rise to a scale-dependent sound speed, that is if 
s 7^ 0, the terms at each order in the expansion also have a scale-dependence^. Explicitly 
writing the scale dependence of the leading term at each order in the Lagrangian 
(Eq. (|3.25|) ) not dependent on A/E and other new parameters, we find 

for n > 3 and where q = {ris — l)/2 — 2s. Different factors of 27r, etc, in the definition 
of will make no difference. Clearly, if s < {ris — l)/4 ~ —0.01, higher order terms 




*See p8, p2, 03] for a discussion of scale-dependence in sfow-roU loop calculations. 
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will have a stronger scale dependence. How can we be sure, then, that the fluctuation 
Lagrangian expanded to order n has actually captured the dominant contributions on 
scales smaller than ko7 This is relevant, for example, for calculations that involve the 
probability density function, which depends in principle on all correlation functions 
and is often treated as an expansion in moments around a Gaussian. We can demand 
that at any scale k, only a finite number of terms in the expansion contribute. Label 
fc* the scale where £„ = 1. Then if /c* > k^_i for all n, we can be sure that at scale fc* 
the observables calculated from the Lagrangian expanded to order n are correct. Then, 
we can write 

/i,*\9("-2) / u* \g{n~3) 

^"1 -'-'{^) '"^ 



For the case n = 3, this gives 



h* 



Then demanding fco/^s < 1 also means that for g > 0, (ko/kl)'^ < 1, and so I/3 < 1. 
But L3 = V]^"^ {kojc'"^ {ko) , so we have just recovered Eq. (|3.10| ). 



C 

For n > 3, we have 



^"l = (4.9) 

, * X q{n-2) / N (n-2)/(n-3) 



Again, the left hand side must be less than one, so finally we have 
The implications of this result are shown graphically in Figure 0. 



5. Multiple Fields 

In this section we will generalize the bounds to the case of multiple fields. This is not 
as precise as we were able to make it in the context of single field slow-roll but the 
basic idea remains the same. This bound would apply to models of assisted inflation 
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Figure 2: (a) Scaling of £3 (solid blue), £4 (long-dashed red), £5 (short-dashed green) with 
coefficients that obey the bound in Eq.( |4.10[ ), with s = —0.1. (b) Same as before, but with 
coefficients that do not obey the bound. Note that when the equality is saturated, all terms 
would meet at the point where C3 crosses 1. 



p9| with scalar fields and the following action (note that in the following the sound 
speed of all the fields are set to 1). 



S = Jd'x^J2{l^^ 



(5.1) 



In the simple case of no cross-terms between the fields, W = '^iV{(f>i). Such models 
were proposed in string theory (see [^, ^ |^ for a recent discussion of 



their viability). Assuming that there is a homogenous and isotropic solution of the 
fields' equations of motion, we can assume an FRW metric with Hubble scale given by 



(5.2) 



Each field whose mass m^^ is smaller than H will undergo quantum perturbations of 
order 50j ^ ~ H where the subscript q stands for quantum. Neglecting perturbations 
of the metric, the action at second order is 



) 



During inflation, the gradient term (and kinetic term) will be of order 



(5.3) 



(5.4) 
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since each field has perturbations of order H over distance scales of order 1/H. De- 
manding that this is smaller than the original background gives the following large N 
bound 

Note that no conditions have been put on any of these scalars fields except that their 
masses are smaller that H so that they fluctuate. So N counts the number of scalar 
fields with low masses in the effective field theory, these fields do not need to have large 
vev or be "part" of the infiaton for this bound to hold. Finally note that this is a very 
simple derivation of a bound that has spurred many recent activities from black hole 
physics [|l^, 13 1 . A similar bound was derived in a different way in [46, 47 1 while |4S] 
used it to put constraint on eternal inflation in Large N models. 

One can also look at higher order interactions in S^. While we expect all sorts of 
couplings between the various fields, there will be, among others 

S^ = j dtdPxa{t)e^'^ 50, + . . . (5.6) 



which gives 



f-eV^E^'^^ (5.7) 

'-'2 , 



But Sipi is a random variable with mean zero and variance H. So the linear sum over 
50j has a mean of zero and a variance of \/NH. Using the variance to estimate the 
value of the contribution from 5*3 gives the following condition: 



6. Conclusions 

In this paper, we have investigated the conditions under which the usual story of 
a slightly perturbed inflationary background, with the variance of the perturbations 
computed from the quadratic action in fluctuations, is valid. Requiring a perturbative 
calculation gives simple bounds on the Hubble scale, sound speed and the number of 
low mass fields. The bound on the Hubble scale is simply 
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which for a small sound speed can require a significantly subplanckian Hubble scale. A 
stronger bound on the sound speed can be found by looking at the higher order terms 
in the action. We found that 

When this bound is violated, perturbation theory breaks down since terms in S'3, 5*4 and 
higher order all become of the same order as the quadratic terms. This is a signal that 
the C field is becoming a strongly interacting field whose variance cannot be computed 
by perturbative methods. This bound also implies that a period of stochastic eternal 
infiation is always outside of the pertubative regime in models with sound speed less 
than 1 (really 1 — > e — 2s, for the terms in Eq. |2.10| ). 



When this bound on is translated into a bound for the magnitude of non- 
Gaussianity from the three-point function, we find that it is less stringent than the 
current experimental bound from CMB data. Nevertheless, for a scale dependent and 
decreasing sound speed, we find that perturbation theory may break down before one 
achieves the necessary number of e-folds for a successful infiationary era. We found 
that the number of e-folds before entering the non-perturbative regime is 



(6.3) 



Without an understanding of what happens in the non-perturbative regime, this can 
severely constrain models like DBI infiation where the sound speed varies. The scale 
dependent constraints (the bound on Cg during the full range of infiation and some 
constraint on higher order interactions) must therefore be included in any analysis of 
this kind of model. 

These constraints should be of interest in the popular program of reconstructing of 
the infiationary action. Very often, no constraints are imposed on the various terms in 
the action that is being reconstructed. We have shown here that there exist conditions 
that may be especially important for reconstructions involving a general kinetic term. 
It would be interesting further work to impose the consistency of perturbation theory in 
the general reconstruction approach in the literature, especially for models with general 
kinetic terms These constraints should also be imposed to restrict the phase space 



of possible configuration when studying such models with Monte-Carlo simulations |^ 

We have been estimating the size of terms in the action to give these bounds. 
Of course, the formal definition of the breakdown of perturbation theory comes from 
looking at the contributions to correlation functions. Although we expect that the 
bounds we have found are similar to demanding that the 1-loop contribution to the 
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two-point function is always small, it would be interesting to show this in detail for the 
general case. We also looked at a simple multi-field model of inflation and used the 
same logic to arrive at an upper bound on the number of fields. 

N<^ (6.4) 

This is exactly the same bound one finds from black hole physics or from renormal- 
ization of the Planck mass. Interestingly the black hole bound seems to indicate that 
the bounds we have found in this paper cannot be violated even in the strong coupling 
regime. Clearly, a more thorough analysis of the multi-fields case is needed and this is 
work in progress. 

Finally, in this paper we have not discussed what happens when perturbation 
theory fails. Does inflation end, does the system dynamically correct, or is there an 
alternate description of the inflationary background? In the case where the higher order 
interactions are small ^3 < 5*2 but the gradient energy dominates over the background 
'S'2 > So, we expect the universe to become inhomogeneous and stop inflating. It would 
be nice to have an explicit description of this transition. A promising way to answer 
these questions is to look at brane inflation where we can hope to use the string theory 
tools to get some handle on the non-perturbative regime. 
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